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Abstract 

In a simple extension of the minimal supersymmetric standard model, out-of-equilibrium decays 
of TeV scale exotic vector-like squarks may generate the baryon asymmetry of the universe. Baryon 
number and CP violation are present in the superpotential, so this mechanism does not rely on 
CP violation in supersymmetry breaking parameters. We discuss phenomenological constraints on 
the model as well as potential signals for the Large Hadron Collider and electronic dipole moment 
experiments. A variation on the TeV scale model allows the exotic squarks to be the messengers 
of gauge mediated supersymmetry breaking. 
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I. INTRODUCTION 



There is a large literature on TeV-scale exotic quarks and squarks (see |l| and the 
references therein). In addition to being of immediate interest for collider searches, they 
are present in many supersymmetric (SUSY) grand unified theories (GUT) and in string 
based constructions. The messenger sector of gauge mediated models of supersymmetry 
breaking (GMSB) also contains such matter [J 0, 0]. In this paper, we investigate the 
possibility that decays of vector-like quarks might have sourced the baryon asymmetry of 
the universe (BAU). In contrast to many models of high-scale baryogenesis, this model 
presents the tantalizing possibility of probing the physics of baryon number generation at 
future colliders. 

To the Minimal Supersymmetric Standard Model (MSSM), we add exotic vector-like 
quark supermultiplets with renormalizable baryon number violating interactions in the 
superpotential. The decays of the lightest exotic particle (LXP) source the BAU. For 
concreteness, we focus on scenarios where the LXP is a squark, rather than its fermionic 
partner. While it is possible to generate the BAU with generic TeV-scale masses for the exotic 
squarks, there are also interesting regions of parameter space where the BAU is achieved via 
a resonant enhancement which requires highly degenerate exotic squark masses. 

The CP violation responsible for the BAU is independent of supersymmetry breaking. So 
far the only non-zero phase observed in nature is the Cabibbo-Kobayashi-Maskawa (CKM) 
phase which resides in the superpotential. Constraints on new CP violating processes 
induced by SUSY breaking are stringent. Perhaps nature has only chosen to have large 
CP violating phases in supersymmetric terms 1 . If nature chooses this path, models that rely 
on SUSY breaking phases cannot generate the baryon asymmetry, and a model such as the 
one presented here would be required. 

Previous work has noted the possibility that out-of-equilibrium, baryon number violating, 
superpartner decays might generate the BAU B0.0J , §]. In particular, Dimopoulos and 
Hall [5| used the baryon number violating operator u c d c d c to produce the BAU from the 
decay of MSSM squarks. In contrast to our approach, CP violation in that model derives 
from SUSY breaking terms. 

There exist a variety of phenomenological constraints and future tests of our model. 
While our approach is largely insulated from the phenomenological difficulties associated with 
EDMs, it predicts values which could be seen in future experiments. Unitarity of the CKM 
matrix and D° — D mixing yield important restrictions. For favorable parameters, one could 
perhaps observe baryon number violation directly at the LHC. Cosmological considerations 
require the reheat temperature to be at most O(10 GeV) for the TeV scale exotics, high 



1 We thank M. Luty for discussions of this point. 
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ion , but low enough to 



enough to allow for a thermal weak scale dark matter candidate 
avoid the gravitino and moduli 12, 13, 14] problems. 

In the next section we describe the details of the model. In Sec. II III we calculate the 
asymmetry from the LXP squark decay. In Sec. [IV] we outline the cosmology and related 
constraints. In Sec. [V] we discuss the variation when the exotics are the messengers of 
gauge mediation. In Sec. ED we discuss the low energy observables and collider signatures. 
The appendices discuss model building challenges for degenerate exotic squarks and provide 
explicit estimates of the cosmological rates for Sec. [TV] 



II. THE MODEL 



The relevant matter is the three generations of colored MSSM chiral superfields (u?, d°, 
qi), i = 1 . . .3, supplemented by N families of exotic vector-like quark superfields (D^Di), 
i = 1 . . . N. We concentrate on the model with N = 2, which is the simplest case where this 
mechanism is viable. There is an approximate Z2 "exotic-parity" under which the D and 
D are odd while all other superfields are even. If this parity were exact, the LXP would be 
stable. The decays of the LXP generate the BAU. The superpotential is 

W = WmSSM + W E xotic, (1) 

with 

W Ex otic = g' ijk u\ D s D k + dl ~Dj + C^-j H d qi Dj + M tj A D„ (2) 

where H d is the MSSM Higgs which couples to down-type fields; v d = (H d ); g' is a B- violating 
coupling between the MSSM and exotic sectors; fi' L and fj,' R are (small) exotic-parity violating 
couplings, and M is the mass matrix for the exotics. This superpotential is in a basis where 
Standard Model Yukawa couplings have been diagonalized, and there is no mixing between 
the MSSM and exotic (s)quarks at gauge boson/gaugino vertices. In estimates below we 
use a common exotic-parity violating coupling, // = fi' L = fi' R . Dependence on either n' L 
or n' R follows from Eqs. (jSJ), (00), and (jSJ). Since // <C M for all viable models, we work to 
lowest order in //. We have omitted similar couplings of the form u c d c d c and u c d c D. Such 
5-violating couplings are present after a rotation to eliminate the // terms. The assumption 
that such couplings are negligible prior to rotation motivates the hierarchy in Eqs. (J4]) and 



The superpotential in Eq. (|2J) satisfies two of the Sakharov conditions [15]: B and CP 
violation. In the presence of both g' and //, it is impossible to consistently assign baryon 



2 A Z4 extension of R-parity (which would also be broken by //) can be constructed to impose the vanishing 
of the u c d c d c and u c d c D couplings. 
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number, and there are physical CP violating phases for N > 2. Using field redefinitions 
of the u\ which leave the mass matrix diagonal, one can always make g[- k real for N = 2. 
Phases remain in the \J matrices. As we discuss in Sec. UV[ the out-of-equilibrium condition 
is dictated by the cosmology: we imagine that the late decay of a modulus reheats the 
universe and (over)populates the LXP. 

We demonstrate the diverse phenomenology of the (g', //) parameter space by presenting 
three scenarios which (see Table [I]) 

i. have generic TeV scale masses for the exotic squarks, 

ii. maximize the reheat temperature of the universe, thereby requiring degenerate TeV 
scale squarks, 

iii. identify the exotics with the messengers of GMSB. 

Also shown in Table d are the exotic squark mass and the splitting between the two lightest 
squarks. For the degenerate D mass benchmarks (ii. and iii. above), all splittings are at the 
sub-percent level, which leads to a resonant enhancement of the BAU. 

We assume no CP violation in the SUSY breaking sector, consistent with our philosophy 
that all CP violation comes from superpotential couplings. Additionally,this both simplifies 
the analysis and highlights differences between our model and that of [5j, where the phase 
arises from soft-terms. While we are agnostic about the origin of the M and fi f terms, for the 



large splittings and high reheat parameters in Table [H the Giudice-Masiero mechanism 16 
might be responsible for their origin (perhaps with a loop-factor generating the hierarchy 
between them). 

A SUSY breaking term bmDD splits the squark masses, resulting in light and heavy 
mass eigenstates, (De)i and (Dh)i, where i — 1,2 for N = 2. Unless the non-holomorphic 
contributions to the squark mass are large, there is an exotic squark lighter then the 
exotic quarks. Should the non-holomorphic SUSY-breaking make the LXP a fermion, the 
generation of the asymmetry proceeds in a nearly identical fashion through the decay of the 
exotic quarks. For unity of discussion, we will assume a squark LXP for all benchmarks. 

The mass difference AM 2 = (Mi) 2 . — (Me) 2 between the two lightest exotics has a large 
impact on the size of the BAU gene rated. A near degeneracy yields a resonant enhancement 



of the baryon asymmetry 13, ll8|, ll9j]. We discuss this possibility in detail in Sec. IIHI 
When required, to achieve a degenerate spectrum we assume that some symmetry enforces 
degenerate values of M for the two families of exotics. The symmetry can be broken by the 
// terms, so we assume an 0(1) generation dependence in fjf^. This implies that the lightest 
(heaviest) two exotic squarks, (De)i ((Dh)i) are degenerate in mass up to jj! 2 corrections. 
This motivates the parametrization of AM 2 in Table [I] Since each is an equal admixture 
of D and D , there are additional factors of 1/ \[2 introduced into the interactions relevant 
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Large D mass splittings 


scenario 


g' 


M 
(GeV) 


(GeV) 


AM 2 
(GeV) 2 


e 


large splittings 


0.4 


500 


4 


(100 /i') 2 


2 x 10~ 6 


Degenerate D masses 


scenario 


g 1 


M 
(GeV) 


(GeV) 


AM 2 
(GeV) 2 


e 


high T RH 
gauge mediation 


0.005 
0.01 


1000 
10 6 


2 
1 


(1.3 M ') 2 

(i.o M ') 2 


8 x 10" 6 
6 x 10~ 5 



TABLE I: Benchmark parameters for three different scenarios. The first has generic values for the 
lightest two D squarks. The last two rely on degenerate squarks to give a resonant enhancement of e: 
a high reheat temperature scenario, and a case where the exotics are identified with the messengers 
of gauge mediation. In the last two cases, we assume all the a'-, and u'.- are approximately 
independent of family. For the first (large splittings) benchmark we require hierarchies in these 
values to avoid the phenomenological bounds alluded to in Sec. Ill Bt 

— * 

for the exotic decays. For clarity, we work with the "helicity" squark eigenstates, D and D . 
We refer to both (De)s as LXPs since both can potentially contribute to the BAU. 

The degeneracy between the two LXP states can be broken by either off-diagonal elements 
in the SUSY and SUSY-breaking masses or by radiative corrections. Since degeneracy is 
important for benchmark points ii. and iii., it is important that these terms can be made 
small. The absence of these terms can be understood in terms of the same (almost) conserved 
family symmetry mentioned above. Some relevant model building issues are discussed in 
Appendix [A] To keep expressions simple we will often use M for both the SUSY mass 
parameter and the mass of the LXP when estimating various processes. 



A. Interactions Induced by Diagonalization 

Mass diagonalization mixes d^ with D and d,R with D. We do not introduce distinct 
notation for gauge and mass eigenstate fields. In what follows, cZl, for example, is the 
mass eigenstate state with the largest overlap with the from above. Because di and D 
have different electroweak charges, the rotation to the mass eigenbasis induces the following 
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off-diagonal couplings to the Z° boson 

where g w is the SU(2) coupling constant and c w = cos9 w . There are also MSSM-exotic 
couplings with the W , as well as the supersymmetric analogues of both of these interactions. 

There is one other class of interactions important for this study. After the d^ — D rotation, 
the following couplings appear: 

9' Mf) u c d c D (4) 

g' {^ju c d c d c . (5) 

The generation of the BAU relies on the interactions in Eqs. ([3]) and (j3J). If either fi' L or /j,' R 
were to vanish, then the BAU generated would be suppressed by powers of ma/M, where rrid 
is a Standard Model down-type quark mass. The interaction of Eq. (jSJ) is a standard MSSM 
R-parity violating coupling. It is small because fi'/M <C 1. For the purposes of this model, 
the dominant effect of these operators is to cause the LSP to decay, which is relevant both 
for collider signatures and for the cosmology. 

We assume off-diagonal gauge interactions between the exotic and MSSM sectors induced 
by soft-terms are sub-dominant for the purposes of the calculation of the BAU. We assume 
this both for simplicity and because we wish to emphasize that this mechanism can occur 
independent of SUSY breaking. 



B. Textures in g' and // 

One possibility is that all g'^ k (and all the /i^-) are comparable. Under this assumption, 
phenomenological constraints make it difficult to realize the BAU without a resonant 
enhancement due to degenerate squarks (more on this in Sec. IHIj) . However, these constraints 
(e.g. D° — D mixing) only pertain to specific families, and can be avoided if hierarchies 
exist in these couplings. We assume such textures apply for the "large splittings scenario" 
of Table [H In this case, the values shown in Table [J are the biggest entries in the g' and // 
matrices - they lead to the dominant contribution to the BAU. 

There is another potential motivation for textures in //. The lightest eigenvalues of the 
full quark mass matrix, which would correspond to the MSSM down-type quarks, are given 
by a see-saw: 

physical = md _^R_ (6) 
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Without textures in the // matrices, M = 500 GeV implies that y! < 2 GeV to avoid fine- 
tuning between the two contributions to the down quark mass. Hence, for the large splittings 
benchmark we assume there is a texture which eliminates this tuning for the down quark. 
This can be done without eliminating all of the CP violating phases. For the benchmarks 
with degenerate LXPs \J is already small enough to avoid fine-tuning. 



III. BARYOGENESIS FROM EXOTIC SQUARK DECAYS 



Our goal is to reproduce the BAU, accurately measured by the WMAP5 [20J data to be 



77 = UB Ub = 6.225 ± 0.170 x 10- 10 , (7) 



where ng (n^) is the number density of baryons (anti-baryons) and n 7 is the number density 
of photons in the universe. 

A calculation of the baryon asymmetry requires two ingredients: a knowledge of the co- 
moving number density of exotics n^/s when they decay out-of-equilibrium and the baryon 
asymmetry created in each decay (e). A necessary condition for the squarks to be out-of- 
equilibrium is that the annihilations are no longer effective. When annihilations eventually 
do freeze-out, the resulting riu/s produces an insufficient rj. Hence, the squarks must be 
populated by some non-thermal source. For our benchmarks, non-thermal decays of a heavy 
field (see Sec. HV|) generate an n D /s in the range 10~ 4 — 10~ 6 , necessitating an e ~ 10~ 4 — 10~ 6 
to reproduce the measured rj. In this section, we discuss the calculation of e and postpone a 
detailed discussion of n D /s to the next section. 



The LXP decays yield a net baryon number per decay 

T(( D e )i ->u + d)- T((D* e )j = m^+ g) 

I r to tal((-D^)i) 



E 



where n is the number of squarks that make non-trivial contributions. Since I?-violating 
decays to the MSSM states are suppressed by // 2 , exotic states have small partial widths for 
these processes. However, exotic-parity ensures the total width of the lightest exotics are 
also suppressed by yf 2 . This allows the B- violating decays for LXP states to compete with 
the total width, yielding an e of appreciable size. 

First we estimate e neglecting possible complications due to resonance. We consider the 
decays in Fig. [Tp. Including the sum over all possible quark final states gives a width: 

rtotai 9£' 2 +6(w(2eg) 2 1 /y^"- 



16 vr 2 M 



~ T7 M. (9) 



3 There are a number of sources of order one uncertainty in our estimates. Since we are not interested in a 
detailed exploration of the parameter space, but rather a demonstration of the viability of this approach, 



7 



D 




+ 



D 




+ 



D 



/ d 




+ 



FIG. 1: Representative contributions to r» otal for the lightest D squarks. For the purposes of our 



estimates we assume the third process is kinematically allowed for all three families of d. 




FIG. 2: Representative diagrams that interfere to give a net baryon number. The squark mass 
insertions illustrate which "helicity" component of the (De)i is contributing to the amplitude. 

A representative contribution to e is given in Fig. [2J We assume D can decay to (three 
families of) d and Z° as well, so there will also be intermediate loops involving these particles, 
which we include in our estimate of e. We estimate 



2 \ Cm J 



-no— resonance 



(10) 



where we have included contributions from both of the light Ds, and defined AM 2 = {Mi)\ — 
{Mi)\. In this estimate, we assume all // and g n s are the same order, and assume a single 
non-zero (maximal) phase. Motivated by phenomenological constraints discussed below, we 
take /i'/M ~ 10~ 2 , g' ~ 1(T 2 which gives 



///mV / m 2 
lo- 2 y v io- 2 J 1am 2 



- no — resonance 



we will not worry about errors of this size - a small change in the input parameters can compensate. In 
this spirit, we display (and compute) representative contributions to e, but do not make an exhaustive 
calculation of all diagrams. For example, we neglect processes involving the and their superpartners 
and LXP decays to u + d. We also assume that Z° is a mass eigenstate and neglect corrections due to 
neutralino mixing. 
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Since we are trying to achieve an e > 10~ 6 , it is clear why a resonant enhancement is 
necessary for much of the parameter space. The need for resonance can only be avoided for 
larger g', which requires a texture to avoid flavor changing neutral current constraints. 



A. Full Calculation of e Including Resonance Effect 

We now calculate e accounting for the possibility of highly degenerate (De)s following 



181 Il9l | . Using the couplings of Eqs. d2J) and (|3D gives 



fe) 2 M(</%i G?')ii2 W k2 Whi W P 2 04 



X 



(m 2 - m|) rr 1 f 



(M 2 -M 2 ) 2 +(rr al § 



-(1<-2U, (12) 



where Mj = (Mt)i, Mi = Ma and a sum over j, fc, p = 1 . . . 3 is implied. We assume the 
physical phase in one of the superpotential couplings (g' A^ijfcC 1 ^ dj D^) equals n/2. The 
value of e for each benchmark is given in Table [H 

We parametrize the mass splitting by (x//) 2 = {{Mi) 2 , — {Mf) 2 ) <C M 2 . Then the \J 
dependence essentially cancels. This illustrates how the resonance effect can compensate for 
small values of fi' . 

We have plotted e in Fig. [3] for the high reheat parameters. Note that e — ► as x — > 0: 
in the limit that the squark masses are degenerate, (De)i is indistinguishable from (Dfa, 
and it is not meaningful to interfere the two. It reaches a maximum e max ~ 9 x 10~ 5 for 
x = 0.1. This plot has been made using a relatively small value of g' = 5 x 10~ 3 . Larger 
values of g' would allow even larger values of e max . Therefore, if splittings are of order /i' 2 , i.e. 
x ~ 0(1), we can achieve e ^> 10~ 10 , as required to generate the baryon asymmetry. One 
also has to keep track of which Dg states contribute to e. In Eq. ( JT2|) we have assumed that 
both r^ otal and rlf* 211 are approximately given by r£ 0tal in Eq. (Q which implies both (Dt)i 
and (De)2 decays can contribute. However, this is only valid when the following decays are 
kinematically forbidden: {pi) 2 ~^ D\+u and (De)2 — > {Dt)\ + x° + u, where x° is the lightest 
neutralino. We assume the first decay is forbidden, as determined by the relative size of M 
and M. For the high reheat parameters the second decay opens up at x m 230 for m^o = 100 
GeV. Once this decay channel opens, the contribution to e from {Di)2 becomes negligible 
and the (1 <-> 2) portion of the curly braces in Eq. (fT2l) drops out. This accounts for the 
apparent discontinuity in Fig. [3] at x = 230; this is the point when the (De)2 contributions 
no longer contribute. Then e ~ 1/x 2 as x — > 00. 
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FIG. 3: Plot of e as a function of the D mass splitting parameter x for the high reheat benchmark 
(see Table U). Note that 

e max ~ 9 x 10 ^ for x — 0.1. We have also taken into account the change 
in r2° tal at x = 230 when new on-shell decay channels open up (see text for discussion). 

IV. COSMOLOGY 

If the universe followed a thermal history, for M < 10 13 GeV, annihilations would keep the 
squarks in equilibrium until they became non-relativistic. Assuming they lived long enough 
to satisfy the out-of-equilibrium condition, the squarks ultimately would have frozen-out with 
ri£>/s <C 7], so their decays could not have generated the BAU. Therefore, there must be a 
non-thermal source for the LXPs. The relevant cosmology begins with a universe dominated 
by a long-lived state, <f), which could be the irmaton or some other modulus. Its decays 
populate the LXPs. 



A. Asymmetry Generated in Decay 

Given a branching ratio (BR) for <fi decaying into the exotic squarks, we approximate the 
co-moving abundance of the LXPs as [2l| : 




where rip is the number density of LXPs; s is the entropy density of the universe; Trh is the 
temperature of the universe generated by the decays, and is the mass of the field. 
Once the exotics are produced, they must decay before they annihilate back to equilibrium, 





(13) 
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i.e., 

P decay > r ann (T RH ), (14) 

where the annihilation rate at the reheat temperature r ann (TRH) is 

2 7T 2 T 4 

r a nn(T RH ) = Y D s (a anQ v) = BRg* — {a ann v). (15) 

45 rritf, 

Here, g* is the number of relativistic degrees of freedom. For this model, one important 
annihilation process is D + D* — > g + g where g is a gluon. We estimate the thermally 
averaged annihilation cross section as 

<'-"> M (£)' (16) 

where g s is the strong force coupling constant. 
We also check 

J- decay 

> H(T RB ), (17) 

where if (T R h) is the Hubble rate evaluated at the reheat temperature. This means the 
exotic squarks decay "instantaneously," i.e. when the temperature is still T RH - 

Then the generated baryon asymmetry is given by 

s = (nn\ = 7MeBR /7rh\ 

ra 7 \ s / \ / 

where e parametrizes the amount of baryon number violation produced by each exotic decay 
(see Sec. [HHfor the calculation of e) and the factor of 7.04 is from the ratio (s/n 7 ) to da y - The 
last task is to make sure that this BAU survives to the present day. 



B. Washout Processes 



There are processes which can washout the BAU. Requiring them to be ineffective 
constrains the maximum reheat temperature. Examples of the most dangerous of these 
baryon number violating processes are 

u + d -> D, (19) 
u + Z° <-> <$ + d\ (20) 
u + g «-> D + D, (21) 

where we assume that Z° is the LSP. The dominant washout process involving only MSSM 
states (Eq. ([2"0~]) ) will always include the LSP, since the LSP suffers the least Boltzmann 
suppression. The process in Eq. (fill is known as inverse decay (ID) and is proportional to 
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Exp(— M/Trh)- The rate for the process in Eq. (12TJ1) receives the Boltzmann suppression for 
the one heavy initial state (Exp(— rhz° /T-bm)) when the process goes from left to right. Due 
to mixing effects (see Eq. (j5J)) the cross section gets additional suppression (cr ~ (/i'/M) 4 ). 
Since the final states are effectively massless at T RH , we will refer to these processes as "light." 
For the GMSB benchmark, T R h 3> ^lsp ; so there will be no Boltzmann suppression. In 
this case the rate will be negligible due to (/z'/M) 4 suppression (see Appendix[B]for details). 
The thermally averaged cross sections for diagrams like those in Eq. (1211) do not depend on 
// but do suffer Boltzmann suppression for T R h < M due to the heavy final states when the 
process goes from left to right. Therefore, we will refer to these processes as "heavy." 



Z° 'X/XyX/XJ 



D 



9 ^-Q~6~6^S~ < 



D 



D 



FIG. 4: Example diagrams for washout processes. 

In Appendix [B] we can estimate these rates. No washout of the baryon asymmetry occurs 
as long as 4 : 

H(T RU ) > r ID (T RH ), (22) 
H(T R u) > rtr h y out (T RH ), (23) 
#(7rh) > rS out (T RH ). (24) 

Ensuring these inequalities places constraints on the four-dimensional parameter space 
spanned by g', fi', T R h, and m^. For various values of the parameters each of the different 
cosmological constraints of Eqs. (fl4l) . (1221) . (123!) and (124")) can become the most important. 
As an illustration, we have shown the allowed region for the reheat temperature as a function 
of the LXP mass for different /x's in Fig. [51 While not the case for the g' chosen in Fig. [51 
for larger g', r^^ ut can be the strongest constraint. 

Once we have constrained T RH and /i', we need to set appropriate values for and BR, 
which enter the expression for Yd, Eq. (fl3|) . One constraint on the inflaton mass is > 2 M 
so that the decays described in Sec. IIVAI will be kinematically allowed. Then the need for 
Y D ~ T RH /m0 to be large competes with keeping r ann ~ T^/m^ small enough to satisfy Eq. 
(JHJ). We chose BR = 0.1 as a reasonable estimate 5 for the branching ratio of — > D^Di . 



4 If these constraints hold at T — Trh then they will hold for all subsequent temperatures since the 
Boltzmann suppression will always dominate over the T 2 dependence of H . 

5 If (f> is a dilaton (i.e. if it enters the Kahler potential in the "no-scale" form) then it will couple to other 
fields in the model proportional to their mass. This can motivate large BRs. 
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FIG. 5: Plots of maximum allowed value for Trh as a function of M for fj,' = 0.1, 1, 10 GeV. The 
shaded region is allowed. The curves bound the regions where H(Tjm) > rJ^^ ut (Tjm) (solid), 
H(T RB ) > r^ out (T RH ) (dotted), H(T RR ) > r ID (T RH ) (dashed) and r decay > r ann (T RH ) (dash- 
dotted). Different processes constrain the reheat temperature for the various values of //. We have 
taken g' = 0.07 and m^ = 2x 10 4 GeV for the purpose of illustration. 

Note that both Eqs. ffTBT) and ffTBT) depend on BR/m^,, so a smaller branching ratio can be 
offset by a smaller inflaton mass. 

Cosmological parameters for the three benchmark scenarios are shown in Table HT1 Exotic 
squark decays can generate the BAU for a wide range of parameters. 



V. GAUGE MEDIATION 



In the simplest gauge mediated models the messengers form N complete 5 and 5 
representations of SU(5). This includes iV families of new vector-like down- type quarks, 
the matter content of our model. Naively, GMSB models exhibit an exact "messenger- 
parity," akin to our "exotic-parity," which in principle could lead to undesirable long-lived 
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Large D mass splittings 


scenario 


GeV 


H 

GeV 


Tdccay 

GeV 


Tann 

GeV 


r m 

GeV 


phcavy 
washout 

GeV 


plight 

washout 

GcV 


GeV 


large splittings 


18 


4 X 10~ 16 


8 x 10" 4 


3 x 10~ 6 


8 x 10" 17 


2 x 10- 24 


1 x 10~ 17 


5000 


Degenerate D masses 


scenario 


Trh 
GeV 


H 

GeV 


1 - decay 

GeV 


Tann 

GeV 


r m 

GeV 


phcavy 
washout 

GeV 


plight 

washout 

GcV 


GeV 


high T RH 
gauge mediation 


75 

1000 


8 X 10~ 15 

1 x 10~ 12 


4 x 10~ 5 
1 x 10~ 8 


1 X 10~ 5 

1 x 10~ 9 


1 X 10~ 15 
~ 


4 x 10~ 16 
~ 


6 x 10" 20 
~ 


10 5 
10 7 



TABLE II: The numerical values for cosmological constraints described in Sec. IIVI for the 
benchmark parameters of Table HI All rates are evaluated at Trh- For the given parameters, 
the numerical value for rj matches the WMAP5 measurement within the approximations made 
here. We expect an 0(1) uncertainty due to diagrams that we have neglected. We have taken 
BR = 0.1, 5*(T RH ) = 100, and m z o = 100 GeV. 



relics (references |22|, |23|, that address this issue). One motivation for models of the type 
presented here is to address this messenger cosmology while simultaneously generating the 
baryon asymmetry. In our models, the low reheat temperature avoids the thermal production 
of messengers. They are instead produced in the decays of the (ft field, and they subsequently 
undergo baryon number violating decays. 

A complete discussion of the messenger cosmology would also require a history for L- 
type messengers, about which we remain agnostic. The lightest one could be a dark matter 



candidate [22| or perhaps there are additional couplings which allow them to decay before 
the onset of big bang nucleosynthesis (BBN), see Sec. IVIDi Either way we assume they do 
not affect the BAU. The final set of benchmark parameters in Tables [J and [TT] is appropriate 
for GMSB. 

To implement the GMSB scenario, we replace the mass term for the exotics in the 
superpotential with 

WGMSB^XSyDiDj, (25) 

where X is a spurion that gets a SUSY breaking vev, (X) = M x + 2 F X - We assume 
identical couplings of different generations to the X field to ensure degeneracy at this order. 
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A. New Contributions to Up-squark Masses 



The tree-level interaction between the messengers and the MSSM u c fields via the 
g' coupling induces new contributions to the up-squark masses from D loops. These 
contributions could potentially spoil the flavor-diagonal nature of the gauge mediated 
couplings. Typically, the leading contribution is at two-loops and is 



At 1-loop, there is an accidental cancellation at O(FxfMx) 2 analogous to |25|, |26j. The 
residual contribution is 

1 F^ 

(&m UR oop ) ij = -— - g ikm g* km (27) 

For the GMSB parameters of Table U 5m 2 - loop » -65 GeV and dm 1 ^ 00 ? » -7 GeV, where 
we have taken Fx = 7 x 10 10 GeV 2 which implies msusY ~ 600 GeV. In the language of 271 ]. 
this leads to a flavor off-diagonal mass-insertions of size 

s»~ . (y "T> 2 - £, (28) 

rl — loop _ (^ 1 - 1 ° OP ) 2 „ , R * 9' 2 Fl (c)Q s 

o = — 2 ^16tt— — r , (29) 

m 9sm M x 

where #sm is the appropriate SM coupling constant. Since these flavor violating contributions 
are in the up sector, the strongest constraint is 5^ O(10~ 2 ) due to D° — D mixing 28]. For 
the GMSB benchmark parameter choices, this provides a (mild) constraint on g', independent 
of \il . 



B. Proton Decay 



It is often stated that both baryon number violation and lepton number violation are 
necessary for proton decay. This is true only if there are no non-leptonic fermions lighter 
than the proton. In GMSB where a light gravitino is present, the decay p — > G + K + is 
open. Following 29j one can estimate the lifetime of the proton in these models to constrain 
the parameters: 

^<0.3GeV(^— j (^eoOGevJ CleV>) v Tb^GevJ- (30) 

We have imposed the proton lifetime to be greater than 2 x 10 33 years [30j for this channel 
using the bound on p — > K + u. For the GMSB parameters in Table [U proton decay constrains 
// ^ 1.2 GeV where we have taken uiq = 16 eV and m SU sY = 600 GeV, corresponding to 
Fx = 7 x 10 10 GeV 2 . Future experiments could discover proton decay if this model is correct. 
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VI. PHENOMENOLOGY 



We begin this section by discussing new contributions to an assortment of low energy 
processes. The need to avoid large violations of CKM unitarity will restrict fi'/M. We 
will find that charmed meson mixing constrains the allowed values of g'. Contributions to 
electric dipole moments could allow a measurement of the // phases in upcoming experiments. 
Finally, we will outline potential collider observables. Neutron-anti-neutron oscillation 
bounds are not relevant in this model due to // suppression. 



A. New Contributions to the CKM Matrix 



The diagonalization procedure of Sec. Ill Al also introduces interactions with the itj, Di 
and W^. This leads to a 3 x (3 + N) CKM-like matrix where the interactions with the 
exotic squarks are suppressed by fi' /M. Using unitarity measurements on the Standard 
Model sector of this new matrix, one can constrain the allowed values of p! jM. The most 



constraining measurement comes from 31 



\VJC + \VusV + \V ub \ 2 = 0.9999 ± 0.001. (31) 

In our model there is an additional contribution of N (/i'/M) 2 on the left-hand side of this 
equation. Requiring this to be within 2 a of the measurement implies 

^ < 0.03, (32) 

where we have assumed N = 2. 



B. Flavor Changing Neutral Currents 



Tree level interactions between the u quarks and the exotic sector with no // suppression 
(Eq. ([2])) give a potentially large contribution to D° meson mixing. The two dominant 
diagrams which contribute to AM D o are shown in Fig. Following [28] this translates into 
a constraint on g'. The experimental limit is Xpo < 8.7 x 10~ 3 with 



AM 



X D o = 



D° 



3.0 x 10 



_ 3 g' A B D o f^ m D o 



ni 



do 



«) 4 + 2(f) 2 Lo, 



D° 



M 2 r 



lm, 



+ m u y 



1 1 M 



M 
M 



(33) 
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where B D o = 0.82, f D o = 0.223 GeV, m D o is the D° mass; m u is the up quark mass; m c is 
the charm quark mass; T^o is the D° decay width; M is the D mass, and M is the D mass. 
This is the source of the constraint on g' in Tabled] for the second benchmark point. Since 
this process only involves the first two generations, this bound can be mitigated by assuming 
a texture for the g' ijk matrix. 



D 



D 



+ 



D 



D 



D 



D 



FIG. 6: New contributions to D° — D° mixing. 

Contributions to neutral kaon mixing are small. These processes are mu' suppressed since 
they involve d quarks. The dominant contributions are from several box diagrams which are 
all strongly suppressed. For the high reheat benchmark parameters (the most dangerous 
case) the contribution to AM K o is roughly 4 orders_of magnitude below the observed value 



31] . Contributions to e' K from AS = 1 operators [32j which only suffer (/i'/M) 2 suppression 



are satisfied for all the benchmarks in Table [H 



C. Electric Dipole Moments 



In principle, there are contributions to electric dipole moments (EDM) due to the phase 
responsible for the BAU. Naively, an 0(1) phase at the TeV scale is dangerous. However, 
in this model, the contributions to EDMs are suppressed by powers of fi'/M. If the exotics 
are the messengers of GMSB, the smallness of this ratio suppresses EDMs beyond anything 
that would be measured. On the other hand, if the exotics are at the TeV scale, the EDMs 
are rendered small enough to avoid current bounds but could be generated at an interesting 
level. 

To generate a non-vanishing EDM requires two elements beyond CP violation: flavor 
mixing and left-right mixing. There are two types of contributions (see Fig. [7]). We consider 
each in turn. 

First, we examine the left diagram of Fig. [71 The left vertex uses the interaction of 
Eq. To close the diagram requires the analytic SUSY breaking insertion d^D. At 
present, the most stringent constraints on hadronic EDMs come from the bound on the EDM 
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(7, flO (7, 9) 



- X - - X 

D /' (d R )j (d L )i \ {dn)j 

1 \ 1 \ 

(d L )j ► ' ^OO^ 1 « (d R ) s {d L ) 3 ► ' ^^00^^ ' fa). 



FIG. 7: Examples of 1-loop contributions to EDMs. The squark mass insertions are b^i and Sfj R , 



i 7^ j, for the first and second diagrams respectively. 



in mercury [331] , du g < 3.1 x 10 29 e-cm. Following 34]. we estimate 



Hg 

Hg „ 2.6 x 10- e cm (Mf) C-^f-) (™™L) 2 , (34) 

V 10 " 3 / V !0" 2 J\ M J ViOOGeVy' V ; 



Requiring the EDM to be below the current bound effectively places a constraint on b^. For 
the first benchmark point, there are no other requirements on the size of this coupling 6 . So, 
this class of diagrams could contribute to an EDM very close in size to the current bound. We 
have employed the QCD sum rules approach to calculating the EDM. In this case, the relevant 
fi' corresponds to mixing with the down quark. This is constrained by considerations of fine- 
tuning, see Eq. ([H]). For cases where degeneracy is important in generating the BAU, is 
constrained by the requirement that enforces the degeneracy of the squarks (see Appendix 
|A]for a detailed discussion). The EDM arising from this class of diagram will be small when 
a resonant enhancement of e is required. 

There is another diagram that can give contributions of an interesting size. For the second 
diagram in Fig. [?l the vertex on the left comes from an interaction induced by the rotations 
to eliminate \J L described in Sec. Ill Al To see the importance of this rotation, note that 
in the SUSY limit there is a correction to the di — d L — Z° coupling at 0{fi' 2 ). While 
hermiticity of the Lagrangian ensures the reality of diagonal interactions with gauge bosons, 
new off-diagonal {di)i — (rfi)j — Z° interactions can inherit CP violation from fi' L . It is the 
super symmetric analog of this interaction that appears in the left vertex. To couple to d^ 
and close the loop then requires an analytic off-diagonal insertion. This leads to an EDM 
for Mercury of the size 

where 5± R is the flavor off-diagonal mass insertion. While this contribution appears small, 
this may be a artifact of the QCD sum rules approach used in the calculation of the mercury 



6 For ba< ~ n' ' M there are modifications of e comparable to the ones calculated here. 
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EDM. Parton quark model 35| calculations of the neutron EDM indicate that the strange 
quark can yield a large contribution. Assuming that this also holds for the Hg nucleus, 
we can make the replacement 5u — > S23, which is only bounded by measurements of the 
branching ratio for b — > 57, implying a constraint 5^ ^ 10 -2 . Also, a ix' L /M as large as 
10~ 2 for the strange quark is allowed without fine-tuning the strange quark mass. All told, 
a Hg EDM as large as O(10~ 30 ) e-cm might be produced by this class of diagram, which 
might be visible in future experiments. Experiments searching for EDMs of the deuteron or 
neutron might also be sensitive to the phases in //. 

Depending on the size of the supersymmetry breaking parameters, it may be possible to 
observe the CP violation responsible for generating the BAU. However, the BAU does not 
depend on supersymmetry breaking, so it is possible that the EDMs might be unobservable 
without affecting the size of the generated asymmetry. 



D. Dark Matter 



Due to the R-parity violating interactions in Eq. (JSJ), the LSP can decay. For neutralinos 
the decay channel will be x° ~^ udd with a rate 



m d 



where g-^ is the coupling between the neutralino, the down-type quark, and down-type 
squark; m x o is the neutralino mass, and fh d is a down-type squark mass. In the large splittings 
(i) and high reheat (ii) scenarios the lifetimes are 0(1O -11 s) and O(10~ 5 s) respectively - 
short enough to decay before BBN. In general, the requirement that the LSP decay before 
BBN places a weak lower bound on fi'/M. For the gauge mediated benchmark (iii), the 
gravitino is the LSP. In this case, the relevant cosmological constraints are for the NLSP, 



which decays via 1/F suppressed couplings (see [36| for a discussion). Note that even though 
there is R-parity violation, the gravitino is stable since it is lighter than the proton and lepton 
number is conserved. 

Since the LSP decays for benchmarks (i) and (ii), it is no longer a viable dark matter (DM) 
candidate. An additional stable weak scale particle must be present. Note, the reheating 
temperature is sufficiently high that this particle (whatever its identity) may have a thermal 
history. If the exotic quarks are embedded into a 5 and 5 representation of SU (5), the lepton 
doublets (L, L) could fulfill this role. However, the simplest DM models with these quantum 
numbers are ruled out due to coherent scattering off nuclei in direct detection experiments. 



19 



E. Collider Signatures 



There are a number of LHC signatures whose observation would lend support to this 
mechanism as the source of the BAU. After confirming SUSY, one would need to determine 
the existence of the exotic quarks or squarks. Exotic heavy colored states at the LHC without 
baryon number violating interactions have been investigated in [l|, [37I, 38, 39 1 , so we will not 



discuss issues related to discovery any further. Instead, we concentrate on phenomenology 
induced by the baryon-violating couplings. 

Much of the novel phenomenology actually arises in the MSSM sector, and does not rely 
on the direct production of the exotics. In particular, the decay of the LSP opens several 
new possibilities for collider phenomenology. LSPs that would otherwise be constrained 
(squarks and charged leptons) are now possible. All MSSM LSP decays will be suppressed 
by (yu'/M) 4 , and over much of the parameter space can lead to displaced vertices. 

Whether or not the LSP decay gives an observable displaced vertex depends sensitively 
on the identity of the LSP. If MSSM squarks are the LSP, they decay directly to two quarks 
with path length 

0-1 W IP' 3 \ 4 / lTeV \ 



^ 1 mm x \Y)Kum)\-^) ■ (37) 

If the charginos or neutralinos are the LSP, they can decay to 3 jets with path length 

/0.i\ 2 /i-0\V 10- 2 \ 4 / fh q \ 4 /100GeV\ 5 



9' J V<W W/M)J \lTeVj V ™* 7 

Here, g% q q is the neutralino-quark-squark coupling. Slepton and sneutrino LSP decays are 
at least four-body (to conserve lepton number), and such decays are too slow to be observed 
on detector time scales. If instead the gravitino is the LSP, then decays of the NLSP to the 
gravitino are much faster than R-parity violating decays. In this case collider signatures will 
be the same as in ordinary gauge mediation. 

The best way to substantiate this baryogenesis mechanism is to observe baryon number 
violating decays. There are four different possibilities which are favorable for the direct 
observation of i?-violating decays: 



1. 2-body B- violating decays of LXPs, 

2. 2-body B-violating decays of LSP squarks with (or without) displaced vertices, 

3. 3-body 5-violating decays of LSP charginos with long-lived tracks, 

4. 3-body 5-violating decays of LSP neutralinos with displaced vertices. 

To be certain a given decay is 5-violating one must establish the baryon number of 
the parent particle. Alternately, this requires observing two different processes: one which 
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preserves B (or alternately, defines the baryon number of the parent particle) and one which 
violates it. In the case of the LSP (2-4 above), there is only a single type of decay (baryon 
number violating). The baryon number of the LSP must be established by a means other 
than decay. One method is examining other particles present in the cascade down to the 
LSP. A second is by searching for associated production processes of the LSP. For the case 
of exotic squark decays (1), one could imagine measuring both D — > x° + d) and D — > u + d 
which would be definitive proof of 5-violation. 

Finally, we point out decays of the type X — > qq are exotic in their own right. The two 
quark final state is indicative of either baryon violation (as here, where X can either be 
the exotic D or a squark LSP) or the existence of an exotic diquark with baryon number 
two. We now comment on the possibility of observing these decays. Ideally, we would like 
to verify that the final state is quarks, and not a quark-anti-quark pair. This is most easily 
accomplished if the final state consists of a top and bottom quark. The charge of the top 
quark can be determined by the charge of the W ± via its leptonic decays. Then, semi-leptonic 
6-quark decays may be used to determine the charge of the initial 6-quark. Competition 
between semi-leptonic decays of the bottom quark, b — > I + and hadronic decays of the 
6-quark followed by semi-leptonic decays of the charm quark, b — > c + ... — >!+..., dilute 
measurements of the bottom charge. To get some idea of how hard it will be to observe 
.B-violation, we will do a simple estimate of how many events should have the right final 
states to reconstruct the 2-body LXP or LSP squark decays to t + b at the LHC In the 
case of the LSP, it is quite likely that there will be a displaced vertex, which should be 
useful in eliminating any SM background. The true challenges are the combinatorics for 
reconstructing a complicated final state (e.g., tbib), the overall event rates, and the relevant 
tagging efficiencies. 

At y/s = 14 TeV, the squark-squark production cross section (for m g = 1 TeV) is ~ 1 pb 
for rriq = 800 GeV, and ~ 5 pb for rriq = 600 GeV 40J, so we will have O(10 5 ) squark pairs 



produced in 100 fb~ of data. We assume a 100% branching ratio to the S-violating channel 
t + b, which is equivalent to imposing a large hierarchy in the g' couplings. We assume a 40% 
efficiency for 6-tagging. Using the 6-quark semi-leptonic branching fraction of 20% and the 
W ± leptonic branching fraction of 22% (electron and muon inclusive), we find 10 3 baryon 
number violating squark decays will be 6-tagged with the appropriate leptonic final states. 
Whether or not this sample will be sufficient to determine the sign of the decaying particles 
is a detailed experimental question. It does not seem unreasonable that such a measurement 
would be possible. 

Because the presence of a (nearly) degenerate pair of squarks is required for the generation 
of a sufficient asymmetry over much of the parameter space, it would be important to verify 
the presence of this duplicity of states. It appears impossible to resolve the squarks at 
the LHC via simple mass measurements. However, a measurement of the production cross 
section might indicate the presence of the second degenerate squark. 
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VII. CONCLUSIONS AND FUTURE WORK 



The baryon number violating decay of exotic squarks can generate the BAU for a wide 
range of parameters. This scenario is constrained by both cosmological arguments and 
particle experiments. The observation of exotic squarks at the LHC would be the first step 
towards verification of this scenario. The subsequent observation of baryon number violating 
decays would be a smoking gun. Alternately, even if no exotics are observed, if the MSSM 
spectrum is consistent with gauge mediation, the observation of proton decay would be 
supportive of this mechanism. 

The goal of this work is to argue the viability of a simple idea: exotic squarks generated 
the BAU. Detailed exploration of parameter space is left for future work. In addition, 
different assumptions about the pattern of supersymmetry breaking might induce additional 
contributions to the BAU. Also, for some regions of parameter space, the exotic quarks could 
be the LXP. In this case, the generation of the BAU proceeds in a nearly identical fashion. 

When very small mass splittings are required, one would like to see specific examples of 
the family symmetries discussed in Appendix |A] to realize the squark degeneracy. This may 
lead to correlations between a solution to the SUSY flavor problem and the BAU. 

Due to R-parity violation, the LSP is no longer a viable dark matter candidate. It would 
be interesting to explore models which could extend the exotic sector to include a new dark 
matter state, thereby connecting the BAU and dark matter. 

Should exotic squarks, their baryon number violating decays, and non-vanishing electric 
dipole moments all be observed, we will have gone a long way towards establishing the origin 
of the baryon asymmetry of our universe. 
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APPENDIX A: ENGINEERING SMALL LXP SPLITTINGS 



We assume exotic-parity is a remnant of an unspecified family symmetry which acts on 
the exotic sector SUSY and soft-breaking parameters. The /i' coupling explicitly breaks this 
family symmetry along with exotic-parity. In the absence of the //, the squarks would be 
exactly degenerate - all splittings are proportional to //. 

Specifically, we assume the symmetry enforces 

Mij = MSij, (A.l) 
Mii = M5 iv (A.2) 
(b M )ij = h M ky (A. 3) 

This will keep new CP violating phases from being generated in the exotic sector. Since // 
is assumed to be the only parameter which breaks the symmetry, all soft-masses which mix 
the exotics with the MSSM will be proportional to this parameter. 



How large a splitting do we expect between the squarks? We first assume Eqs. flA.ljl 



(1A.3j) hold. The introduction of a non-zero // induces AM 2 ~ // 2 . There are additional 
potential sources of splitting. As an example, assume that the term is non-zero. If 
M ~ t^susy the mass splittings are proportional to b^. If instead there is a large hierarchy 
M 3> msusY, e.g. for the GMSB scenario, then the splittings are proportional to 6 2 ,/M 2 . 
Care must be taken to make sure these splittings induced by £y do not upset the resonance 
condition, i.e. AM 2 ~ fi' 2 . This implies £y ~ p! 2 when M ~ msusY, while for larger M 
(GMSB scenario), the weaker condition b^ ~ p! M is required. 

To summarize, in order to achieve degenerate squarks we assume the mass terms obey 
the following properties: 



• The superpotential mass is proportional to the identity: Mij = M 5ij. 

• The non-analytic and analytic exotic squark mass matrices are proportional to the 
identity: = M 5ij and (bM)ij = bM^ij- 

• The mixing between the MSSM and the exotic sectors (cIrD, diD*, etc.) is 
constrained. While the precise size of the mixing depends on the size of M, it will 
at minimum need to be proportional to // M. 

• The MSSM soft-terms are left unspecified. 



Now that we understand the conditions at tree-level we will discuss the radiative stability 
of these requirements. Non-renormalization of the superpotential ensures that we need only 
worry about contributions to wave function renormalization. If the g' coupling breaks the 
family symmetry, loop effects will introduce scalar mass squared splittings of 0{-^-^M 2 ) 
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which can potentially be larger than fi' 2 . We need to make sure that this contribution 
is either sub-dominant or proportional to the identity matrix in flavor space. The latter 
can be done by imposing specific textures on the g' coupling matrix. The form of these 
textures depends on the number of exotic families. The flavor structure of the wave function 
renormalization at 1-loop comes from 

ZS~sLu(iU* ( A - 4 ) 

where we are neglecting contributions 0(/i'/M) since these will not generate splittings larger 
then the tree-level contributions. For N = 2, the couplings have the form 

g'ijk =f ° 9l \ = g[e jk . (A.5) 

\-g'i o 

Then it is easy to see that Eq. (1A.4|) becomes 



gin {g'wY = (a'k) 2 e ji) = - 2 (#fc)%- ( A -6) 

In the case where N = 2, once we achieve splittings of 0(fj,' 2 ) they will be radiatively stable. 
For N > 2, specific textures on the g' couplings will be required for degeneracy to hold at 
loop level. 



APPENDIX B: ESTIMATES OF WASHOUT RATES 



In this appendix we will estimate the rates for the three types of washout processes 
relevant for this model. We will also discuss sphaleron processes which are relevant for 
scenarios where the squarks are the messengers of GMSB. 



1. Inverse Decays 



For inverse decays in Eq. f[T^]) the relevant rate is given by 41 



ID * D 1 decay ' V^A) 

where Y^ q comes from applying the momentum conserving delta function to the distribution 
functions for u and d and integrating over phase space, and r^ e ~^ olatmg is given by the first 
diagram in Fig. [TJ Using the same level of approximation described in Sec. IHII to calculate 
the 5-violating width gives 

45 f M\ V \_ M/Tw /<W 2 



r ID (T RH ) - - " f- e~^r,u . (R2) 

4V2 7r 7 / 2 £f* VTrh/ \32irM ' 
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2. Light Final States 



To determine the rate for the process in Eq. (1201) . we need the number density of the 
Z°s. LSPs are produced in the cascade decays of all exotics and MSSM superpartners which 
were created by the decays of <fi. If either the annihilation rate or decay rate of the LSPs 
is sufficiently fast, the LSPs "instantaneously" re-thermalize and follow a non-relativistic 
equilibrium Boltzman distribution. Otherwise, one should use the n^ from non-thermal 
production. Which of these number densities is appropriate depends on T RH (in the case of 
annihilation) or the size of the coupling in Eq. (j3j) (in the case of decay). For the benchmark 
models, the LSPs do indeed re-thermalize, so the washout rate is given by: 



2 7r / \ An ) \m) (ml VSY + Tl u y ^ 



where n^o is the non-relativistic number density of Zinos, and the label "light" refers to the 
masses of the final states. We have used a simple estimate for the thermally averaged cross 
section. The factor of 12 accounts for the contributions from all MSSM quarks and squarks. 
We have assumed that the bottom quark is "light" for the reheat temperatures considered 
here. For the GMSB benchmark where Trh 3> m^ , the Zino will be relativistic so that 
n^o ~ T^ H . Since (c^Lhout v ) ~ i^'/M) 4 , the large // suppression keeps this rate negligible 
for the GMSB scenario. 



3. Heavy Final States 



For the process in Eq. (|2ip more care is required. The reason is two-fold: there will be 
Boltzmann suppression since only the tails of the u and g distributions have enough energy 
to create the exotic pair, and the initial states are both relativistic at the time of freeze-out 
so Maxwell-Boltzmann statistics do not apply for the normalization. Modifying the results 



of 42J we find: 

/heavy v 
\° washout u l 



2vr 2 T 



poo 

/ - ImDVSK^mds, (B.5) 

Jam 2 

where K\ is the modified Bessel function and the label "heavy" refers to the mass of 
the final states. The factor of K\ assumes Maxwell-Boltzmann statistics which is a good 
approximation since the integral is evaluated from s = AM 2 to oo, which is larger than Trh 
in this study. 
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To obtain a rate we multiply this thermally averaged cross section by the relativistic 
number density for the gluons: 

n g = 1.2 7T 2 (2 x 8)T 3 , (B.6) 

to get 

^ washout = n a (^washout V ) • (B-7) 

In order to evaluate this numerically we approximate the total washout cross section by 
the dominant process which proceeds via a t-channel D quark (see Fig. H]). In Fig. Owe 
show the maximum allowed reheat temperature as a function of the M for all the washout 
rates. 



4. Sphaleron Processes 



For the GMSB benchmark, Trh > n^w which implies that electroweak sphaleron rates 
are unsuppressed. Our mechanism for generating the baryon asymmetry does not generate 
any associated lepton asymmetry. This implies that we have a non-zero value for B but not 
L (where L is the lepton number of the universe). Since sphalerons only violate B + L, all 
they will do is distribute the B generated from the exotic squark decays to both B and L. 
The result is that the sphalerons will reduce rj by approximately a factor of 2. For the exact 



relationship between the initial B — L and final B and L asymmetries, see 43 
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